The matter Lagrangian and the energy-momentum tensor in modified gravity with 
non-minimal coupling between matter and geometry 



O 

(N 

00 



o 
cr 



(N 
> 

m 
in 

o 
o 



X 



T. HarkcQ 

Department of Physics and Center for Theoretical and Computational Physics, 
The University of Hong Kong, Pok Fu Lam Road, Hong Kong, P. R. China 
(Dated: February 8, 2010) 

We show that in modified f(R) type gravity models with non-minimal coupling between mat- 
ter and geometry, both the matter Lagrangian, and the energy-momentum tensor, are completely 
and uniquely determined by the form of the coupling. This result is obtained by using the varia- 
tional formulation for the derivation of the equations of motion in the modified gravity models with 
geometry-matter coupling, and the Newtonian limit for a fluid obeying a barotropic equation of 
state. The corresponding energy-momentum tensor of the matter in modified gravity models with 
non-minimal coupling is more general than the usual general-relativistic energy-momentum tensor 
for perfect fluids, and it contains a supplementary, equation of state dependent term, which could 
be related to the elastic stresses in the body, or to other forms of internal energy. Therefore, the 
extra-force induced by the coupling between matter and geometry never vanishes as a consequence 
of the thermodynamic properties of the system, or for a specific choice of the matter Lagrangian, 
and it is non-zero in the case of a fluid of dust particles. 

PACS numbers: 04.50.+h,04.20.Cv, 95.35.+d 



I. INTRODUCTION 

A very promising way to explain the recent observa- 
tional data [l|, H| on the acceleration of the Universe 
and on dark matter is to assume that at large scales 
the Einstein gravity model of general relativity breaks 
down, and a more general action describes the gravita- 
tional field. Theoretical models in which the standard 
Einstcin-Hilbert action is replaced by an arbitrary func- 
tion of the Ricci scalar R, first proposed in Q , have been 
extensively investigated lately. Cosmic acceleration can 
be explained by f(R) gravity [3], and the conditions of 
viable cosmological models have been derived in Q. In 
the context of the Solar System regime, severe weak field 
constraints seem to rule out most of the models proposed 
so far [H, 0] , although viable models do exist fc!4Tl| . The 
possibility that the galactic dynamic of massive test par- 
ticles can be understood without the need for dark mat- 
ter was also considered in the framework of f(R) gravity 
models [r3 - |l6l . For a review of f(R) generalized gravity 
models see |17j . 

A generalization of the f(R) gravity theories was pro- 
posed in [18| by including in the theory an explicit cou- 
pling of an arbitrary function of the Ricci scalar R with 
the matter Lagrangian density L m . As a result of the 
coupling the motion of the massive particles is non- 
geodesic, and an extra force, orthogonal to the four- 
velocity, arises. The connections with MOND and the 
Pioneer anomaly were also explored. This model was 
extended to the case of the arbitrary couplings in both 
geometry and matter in (ljjj . The implications of the 
non-minimal coupling on the stellar equilibrium were in- 



vestigated in [20| , where constraints on the coupling were 
also obtained. An inequality which expresses a necessary 
and sufficient condition to avoid the Dolgov-Kawasaki in- 
stability for the model was derived in [2l|. The relation 
between the model with geometry-matter coupling and 
ordinary scalar-tensor gravity, or scalar-tensor theories 
which include non-standard couplings between the scalar 
and matter was studied in (27j ■ In the specific case where 
both the action and the coupling are linear in R the ac- 
tion leads to a theory of gravity which includes higher 
order derivatives of the matter fields without introduc- 
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ing more dynamics in the gravity sector 22j . The equiv- 
alence between a scalar theory and the model with the 
non- minimal coupling of the scalar curvature and matter 
was considered in [2j]. This equivalence allows for the 
calculation of the PPN parameters /3 and 7, which may 
lead to a better understanding of the weak-field limit of 
f(R) theories. The equations of motion of test bodies in 
the nonminimal coupling model by means of a multipolc 
method were derived in [24]. The energy conditions and 
the stability of the model under the Dolgov-Kawasaki cri- 
terion were studied in (25[. For a review of modified f(R) 
gravity with geometry- matter coupling see [26| . 

The extra force in the f(R) gravity with non-minimal 
coupling has the intriguing property that it depends on 
both the matter Lagrangian, and the energy-momentum 
tensor. It was pointed out in (2?| that adopting the stan- 
dard Lagrangian density cl„ — P, where P is the pres- 
sure, the extra force vanishes in the case of dust. Differ- 
ent forms for the matter Lagrangian density C m , and the 

resulting extra- force, were considered in [28| . and it was 

(2) 

shown that the more natural form for C m , C m = —p, 
does not imply the vanishing of the extra-force. The 
impact on the classical equivalence between different La- 
grangian descriptions of a perfect fluid was also analyzed. 
Recently, the problem of the two Lagrangian descriptions 
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By assuming that the Lagrangian density L m of the 
matter depends only on the metric tensor components 
g^v, and not on its derivatives, we obtain = L m g^ v — 
2dL m / dg^ LU . By taking into account the explicit form of 
the field equations one obtains for the covariant diver- 
gence of the energy-momentum tensor the equation [l8j 



V% = 2{V"ln[l + A/ 2 (i?)]} 



of a perfect fluid, and its effect on the extra force, was momentum tensor of the matter as [3C 
reconsidered in [29[ . The main conclusion of this study is 
that" It is a fact that by choosing C\ = P there is no ex- 
tra force on a dust fluid, and it is equally undeniable that 
by choosing £2 = —p there will be such a force, which 
may ultimately provide an alternative to dark matter." 
It is also pointed out that " as soon as this fluid is coupled 
explicitly to gravity ...the two Lagrangian densities cease 
to be equivalent" [29| . 

It is the purpose of the present paper to show that in 
the f(R) gravity with non- minimal coupling the matter 
Lagrangian, and the corresponding energy-momentum 
tensor, are not model-independent quantities, but they 
are completely and uniquely determined by the nature of 
the coupling between matter and geometry. This result 
can be obtained by deriving first the equations of motion 
in the modified gravity model from a variational prin- 
ciple, and then considering the Newtonian limit of the 
particle action for a fluid obeying a barotropic equation 
of state. The corresponding energy-momentum tensor 
of the matter is more general than the usual general- 
relativistic energy-momentum tensor for perfect fluids, 
and it contains a supplementary term that may be related 
to the elastic stresses in the body, or to other sources of 
internal energy. The matter Lagrangian can be expressed 
either in terms of the density, or in terms of the pressure, 
and in both representations the physical description of 
the system is equivalent. Therefore the presence (or ab- 
sence) of the extra-force is independent of the specific 
form of the matter Lagrangian, and it never vanishes, 
except the case of (un)physical systems with zero sound 
speed. In particular, in the case of dust particles, the 
extra-force is always non-zero. 

The present paper is organized as follows. The matter 
Lagrangian and the energy- momentum tensor in f(R) 
gravity with non-minimal coupling are derived in Section 
II. We discuss and conclude our results in Section III. In 
the present paper we use the Landau-Lifshitz (3(J sign 
conventions and definitions, and the natural system of 
units with c = 8nG = 1. 







(2) 



dg 
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As a specific example of generalized gravity models 
with linear matter-geometry coupling, we consider the 
case in which the matter, assumed to be a perfect ther- 
modynamic fluid, obeys a barotropic equation of state, 
with the thermodynamic pressure p being a function of 
the rest mass density of the matter (for short: matter 
density) p only, so that p — p (p). In this case, the matter 
Lagrangian density, which in the general case could be a 
function of both density and pressure, L m — L m (p,p), or 
of only one of the thermodynamic parameters, becomes 
an arbitrary function of the density of the matter p only, 
so that L m = L m (p). Then the energy- momentum ten- 
sor of the matter is given by 



dp 



(4) 



where the four- velocity u M = dx^/ds satisfies the con- 
dition g tlu u ll u u = 1. To obtain Eq. (Q} we have im- 
posed the condition of the conservation of the matter 
current, V ' v {pu v ) — 0, and we have used the relation 
Sp = (1/2) p{g^ v — Ufj,Uu) 5g^ u , whose proof is given in 
the Appendix. With the use of the identity u v V — 
d 2 x^/ds 2 + T^ x u"u x , from Eqs. J3J and (gj we obtain the 
equation of motion of a test fluid in the modified gravity 
model with linear coupling between matter and geometry 



II. MATTER AND ENERGY-MOMENTUM 
TENSOR IN MODIFIED GRAVITY WITH 
LINEAR COUPLING BETWEEN GEOMETRY 
AND MATTER 



as 



where 



ds 2 



(5) 



The action for the modified theories of gravity pro- 
posed in takes the following form 



S = 



fi{R) + [l + Xh(R)]L r . 



-9 d 4 x , (1) 



where fi(R) (with i = 1, 2) are arbitrary functions of the 
Ricci scalar R, and L m is the Lagrangian density cor- 
responding to matter. The strength of the interaction 
between f2{R) and the matter Lagrangian is character- 
ized by a coupling constant A. We define the energy- 



/" = -V„ln [l + \f 2 (R)} 



d L m (p) 
dp 



(6) 

The extra- force / M , generated due to the presence of 
the coupling between matter and geometry, is perpen- 
dicular to the four-velocity, f^u^ = 0. The equation 
of motion Eq. ([5]) can be obtained from the variational 
principle 



SS P = 5 J Lpds = 5 J \fQ \j ' g fiv u^u 1/ ds = 0, 



(7) 
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where S p and L p = \/Q^/g llv u> i u v are the action and 
the Lagrangian density for the test particles, respectively, 
and 



Q = [l + Xf 2 (R)} 



dL rn (p) 

dp 



(8) 



To prove this result we start with the Lagrange equa- 
tions corresponding to the action ([7]), 



d_ ( dLp\ _ dLp _ Q 
ds V du x J dx x 



(9) 



Since dL p /du x — y/Qu\ and dL p /dx x — 
(1/2) y/Qg^.xu^W + (1/2) Q y \/Q, a straightforward cal- 
culation gives the equations of motion of the particle as 



+ T» x u v u x + (uV - gT) V„ In y/Q = 0. (10) 



d 2 x^ 
lb? 



By simple identification with the equation of motion of 
the modified gravity model with linear matter-geometry 
coupling, given by Eq. ([5]), we obtain the explicit form of 
y/Q, as given by Eq. ©. When ^/Q — > 1 we reobtain the 
standard general relativistic equation for geodesic mo- 
tion. 

The variational principle ([7]) can be used to study 
the Newtonian limit of the model. In the lim it of 
the weak gravitational fields, ds ss ^/l + 2(f) — iPdt « 
(l + <fi — v 2 /2) dt, where <f) is the Newtonian potential 
and v is the usual tridimensional velocity of the fluid. 
By representing the function \fQ as 



Q= dL„(p) +XMR) dL m (p) 



dp 



dp 



(11) 



in the first order of approximation the equations of mo- 
tion of the fluid can be obtained from the variational 
principle 



S 



dp 

and are given by 



dL m (p) v xf 2 {R)- Lm ^ 



dp 



'2 i 



dt = 0, 



(12) 



a= -V0- V dL ™^ - VU E = a N + d H +d El (13) 
dp 

where a is the total acceleration of the system, o/v = 
— V</> is the Newtonian gravitational acceleration, and 
a E = -VU E = -AV [f 2 (R)dL m (p) /dp] is a supple- 
mentary acceleration induced due to the coupling be- 
tween matter and geometry. As for the term an — 
—V[dL m (p)/dp\, it has to be identified with the hy- 
drodynamic acceleration term in the perfect fluid Euler 
equation, 



a H = -V 



dLm (p) 
dp 



= -V 



"dpdp 
, d P P 



(14) 



where po, an integration constant, plays the role of a 
limiting density. Hence the matter Lagrangian can be 
obtained by a simple integration as 



L m (p) = p[l + U(p)}- [dp, 
•>pa 



(15) 



where IT (p) = dp/ p, and we have normalized an ar- 
bitrary integration constant to one. pq is an integra- 
tion constant, or a limiting pressure. The corresponding 
energy-momentum tensor of the matter is given by 

= {p [1 + $ (p)] +p(p)} u»u» - p (p) <T , (16) 



respectively, where 



—dp 

PO P 



n(p) 



p(p) 
p 



(17) 



and with all the constant terms included in the defini- 
tion of p. By introducing the energy density of the body 
according to the definition e = p [1 + $ (p)], the energy- 
momentum tensor of a test fluid can be written in the 
modified gravity models with geometry-matter coupling 
in a form similar to the standard general relativistic case, 
= [e (p) + p (p)] u»u v - p (p) 

From a physical point of view $ (p) can be inter- 
preted as the elastic (deformation) potential energy of the 
body, and therefore Eq. (|16p corresponds to the energy- 
momentum tensor of a compressible elastic isotropic sys- 
tem. The matter Lagrangian can also be written in the 
simpler form L m (p) = p<& (p). 

If the pressure does not have a thermodynamic or ra- 
diative component one can take po — 0. If the pres- 
sure is a constant background quantity, independent of 
the density, so that p = po, then L m (p) = p, and the 
energy-momentum tensor of the matter takes the form 
corresponding to dust, T^ w — pu^u v . 



III. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have shown that in the 
f(R) gravity with non- minimal coupling the matter La- 
grangian and the energy-momentum tensor can be ob- 
tained uniquely and consistently from the form of the 
coupling between matter and geometry. The coupling 
completely fixes both the matter Lagrangian, and the 
energy-momentum tensor. Since the matter is supposed 
to obey a barotropic equation of state, this result is inde- 
pendent of the concrete representation of the matter La- 
grangian in terms of the thermodynamic quantities. The 
same results are obtained by assuming L rn = L m (p) - due 
to the equation of state p and p are freely interchange- 
able thermodynamic quantities, and the Lagrangians ex- 
pressed in terms of p and p only are completely equiva- 
lent. More general situations, in which the density and 
pressure are functions of the particle number and tem- 
perature, respectively, and the equation of state is given 
in a parametric form, can be analyzed in a similar way. 
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The form of the matter Lagrangian, and the energy- 
momentum tensor, are strongly dependent on the equa- 
tion of state. For example, if the barotropic equation of 
state is linear, p = (7 — 1) p, 7 = constant, 1 < 7 < 2, 
then 



L m (p) = p ll + (7 - 1) 



In 



po 



- 1 



(18) 



and $ (p) = (7 — 1) In (p/p ), respectively. In the 
case of a polytropic equation of state p = Kp 1+1 / n , 
K, n ^constant, we obtain 



L m {p) = p + K 



n + 1 



1 P 



l + l/n 



(19) 



and <f> (p) — Knp 1+1 / n — np(p), respectively, where we 
have taken for simplicity po = po = 0. For a fluid satis- 
fying the ideal gas equation of state p = kspT/p, where 
ks is Boltzmann's constant, T is the temperature, and p 
is the mean molecular weight, we obtain 



L m (p) = p \ 1 + 



In [ — ) - 1 

Po 



+ Po- (20) 



In the case of a physical system satisfying the ideal gas 
equation of state, the extra-acceleration induced by the 
presence of the non-minimal coupling between matter 
and geometry is given by 



a E 



/ 2 (i?)ln 



P_ 



(21) 



and it is proportional to the temperature of the fluid. It is 
also interesting to note that the limiting density and pres- 
sure po and po generate in the energy-momentum tensor 
some extra constant terms, which may be interpreted as 
the dark energy. 



In conclusion, the extra-force induced by the coupling 
between matter and geometry does not vanish for any 
specific choices of the matter Lagrangian. In the case of 
the dust, with p = 0, the extra force is given by 



/" - -V„ In [1 + Xf 2 (R)} - 9 n , 



(22) 



and it is independent on the thermodynamic properties 
of the system, being completely determined by geometry, 
kinematics and coupling. In the limit of small velocities 
and weak gravitational fields, the extra-acceleration of a 
dust fluid is given by 



a E = -W [f 2 {R)\ 



(23) 



The thermodynamic condition for the van- 
ishing of the extra-force is dL m /dg^ v = 
(1/2) (dL m /dp) pig^v — u^u v ) — only. If the matter 
Lagrangian is written as a function of the pressure, 
then dL m /dp = (dL m /dp) (dp/dp), and for all physical 
systems satisfying an equation of state (or, equivalently, 
for all systems with a non-zero sound velocity), the 
extra-force is non-zero. Therefore, the geometry-matter 
coupling is introduced in the generalized gravity mod- 
els [ll] and in a consistent way. The coupling 
determines all the physical properties of the system, 
including the extra-force, the matter Lagrangian, and 
the energy-momentum tensor, respectively. 
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Appendix: Proof of the relation 

8p = \p (gr M „ — UftUv) Sg^ v 

Let's consider a fluid with density p in a system of 
arbitrary coordinates x x . Let's introduce an arbitrary 
coordinate transformation to a new coordinate system 
o/ 1 , so that x a = f a (a A ), a = 0,1,2,3, A = 0,1,2,3. 
In the new coordinate system a A the components of 
r? MV (a A ) are given by rj^ (a A ) = 



the metric tensor 



g a f) (x x ) {df a /da») (dfP/da v ). We define the four- 
velocities of the fluid in the two frames as u a = dx a /ds( x ) 
and Q a — da a /dsi a \, where ds 2 , x 



(.,) 



9af3 



dx a dx^ and 



ds?s = rj a pda a da^, respectively. The fluid obeys the 
equations of mass continuity, which in the two reference 
frames are given by 



1 



d 



-g dx a 



-gpu 



= 0. 



(24) 

In Eqs. (fMf we have taken into account that p is a scalar 
quantity, that is, it is invariant with respect to this coor- 
dinate transformation. Let's now assume that the coor- 
dinate transformation f a takes our current coordinates 
of the particles in the fluid, x A , to the initial coordinates 
of the particles a A (the Lagrange coordinates) . The cor- 
responding Lagrange coordinate system is a comoving 
coordinate system - in it the fluid appears to be at rest. 
In Lagrange coordinates the space-like coordinates a k are 
constants, a k = a k (0), k = 1,2,3. The fluid velocity in 
the comoving coordinate system is 



Q a = 



da° 



ds 



— On 



O) 



9a/3 



if 



df a df 



da° da? 



-1/2 



(25) 



In the Lagrangian coordinates a A the equation of con- 
tinuity becomes 



d 
da 1 * 



9a/3 (f 



df a df p 



-1/2- 



= 0, (26) 



and can be integrated immediately to give 



5^(/ A )So-Srf F ( al ( )' a2 ( )' a3 (°)) 



da da 

where F is an integration function, depending on the 
Lagrange coordinates a k , k = 1,2,3, only. Taking into 
account that y/—rj = \/—g [d(x)/d(a)], where d{x)/d(a) 
is the absolute value of the Jacobian of the coordinate 
transformation x a — > a a , we obtain 

(28) 

In the a a coordinates the velocity four-vector takes the 
form 



ds 



(«) 



6 



Since the Jacobian d(x)/d(a), as well as the function 
F (a 1 (0), a 2 (0), a 3 (0)) , does not contain the quantities 
g M „, the variation with respect to g^ u in Eq. (|28|) gives 



d(x) . , x 

W) {P ^ ] = 2 



But 



8aP da 



= P\ 



d(a)> 



(30) 



(31) 



and therefore 

S (pV^g) = \^-9 Sg^ = ~V-9 {pu^u v ) 5g^ u 

(32) 

Combining this equation with the known relation 
SV=9 = - (1/2) V^99^Sg^ gives 



S P = \p (flV - u^u y ) 8g^. 



(33) 



